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We investigate the isometric transient response of muscle using a quantitative stochastic model 
of the actomyosin cycle based on the swinging lever-arm hypothesis. We first consider a single 
pair of filaments, and show that when values of parameters such as the lever-arm displacement 
and the crossbridge elasticity are chosen to provide effective energy transduction, the T2 curve 
(the tension recovered immediately after a step displacement) displays a region of negative slope. 
If filament compliance and the discrete nature of the binding sites are taken into account, the 
negative slope is diminished, but not eliminated. This implies that there is an instability in the 
dynamics of individual half-sarcomeres. However, when the symmetric nature of whole sarcom- 
eres is taken into account, filament rearrangement becomes important during the transient: as 
tension is recovered, some half-sarcomeres lengthen while others shorten. This leads to a flat T2 
curve, as observed experimentally. In addition, we investigate the isotonic transient response and 
show that for a range of parameter values the model displays damped oscillations, as recently 
observed in experiments on single muscle fibers. We conclude that it is essential to consider the 
collective dynamics of many sarcomeres, rather than the dynamics of a single pair of filaments, 
when interpreting the transient response of muscle. 



I. INTRODUCTION 

The transient response of muscle to a sudden ad- 
justment of its extension, or to an abrupt change in 
load, has been one of the most important sources of 
information about the mechanism of contraction for 
over three decades. Ever s ince the pioneering work of 
iHuxlev and Simmons! experimental data on tran- 

sients has informed theoretic al models of the interac- 
tion between myosin and actin Jchen and Brenner! 1 199r ; 
[ Eisenberg et all Il98fi 1511 Il974t IHuxlev and Simmons , 
Il97li lHuxlevll2000(l . providing a more detailed picture 
than could be obtained f rom the force-velocity relation 
|ffiil [1939; H uxlevl Il957|) alone. The reason is that the 
actomyosin interaction involves several processes which 
occur on different time scales, and these individual com- 
ponents can_be resolved during the transient re s ponse 
ilBremieit Il99ll iBrenner et all Il995t iFord et all Il977t 
119851 Il986t IHuxlev and Simmonsl Il97l[) . The quickest 
process is the elastic deformation of the myosin cross- 
bridges that link the thick and thin filaments. Rapid 
transitions between two or more different bound states 
of the myosin molecule are thought to be the next fastest 
events, while detachment and reattachment of myosin 
heads occur on a slower time scale. 

In an experiment to determine the isometric transient 
response, a muscle fiber is held at both ends to pre- 
vent it from contracting. The muscle is then suddenly 
shortened (or stretched) by a fixed amount, and the ten- 
sion T that it generates is measured. Immediately after 
the imposed change of length, the tension shifts from 
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the isometric value T to a new value, which is termed 
T\. But shortly afterwards (typically within 2 ms), the 
tension adjusts to a new value, termed T%. Subse- 
quently, it gradually reverts to the original isometric 
value To, and the entire transient response is usually com- 
pleted in a fraction of a second. It is generally accepted 
that the initial response T\ corresponds to the mechan- 
ical deformation of cross- bridges and provides a direc t 
measure of their elasticity ijHuxlev and Simmon si I1971IK 
The interpretation of T2 is rather more controversial. 
It is often attributed to force generatio n by the work- 
ing s t roke o f bound myo s in molecules iBrenner et al. , 
1995 | iDukel Il999l ElOfl lEisenberg et all Il98(l Iffil , 
1974|:lHuxlev and Simmonsl . ll97lHHuxlev and Tideswel , 
1996j) and this interpretation has recently gained sup- 
port from X- ray interference techniques applied t o short - 
cning fibers (|Irving et all l2000t iPiazzesi et all l2002c|) . 
But alternative models suggest that the force regener- 
ation might be due, in part, to the r apid binding of 
new myosin h eads to the thin filament ijBrennerl Il99lt 
iHowardl l200l|) , or that it might involve the activation of 
the second myosin head (|Huxlev and Tideswell Il997j) . 

In this article, we wish to address a fundamental prob- 
lem connected with the interpretation of force transients. 
The present theories are all based on the consideration 
of a single pair of filaments, i.e. one filament containing 
myosin molecules, interacting with one actin filament. 
The dynamics of this filament pair is generalized to that 
of a whole muscle fiber by assuming that all filament 
pairs in a fiber behave in exactly the same way. This 
assumption is certainly justified as long as there are no 
static or dynamic instabilities in the system. However, 
the possibility of such instabilities has been known for 
a long time l|HilJ . 11974). Moreover, a stochastic model 
of the actomyosin cycle, based on the swinging lever-arm 
hypothesis, has shown that instabilities do arise when val- 
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ues of parameters such as the lever-arm displacement and 
the crossbridge elasti city are chosen to provide effective 
energy transduction (jDukel 119991 l200fTft . Such instabili- 
ties would give rise to a region of negative slope in the T 2 
curve of a single filament pair. Several reasons have been 
advanced for the absence of any ne gative slope in the ex- 
perim entally determined T 2 curve. iHuxlev and Simmons! 
i|l97lj) argued that the power stroke is sub-divided into 
several small steps, and fixed the step size so that the 
T 2 curve had zero slope for limi tingly small changes o f 
length. In the model proposed bv lEisenberg et afl (l980), 
the flatness of the T 2 curve was explained by a broad dis- 
tribution of cross-bridge strain after attachment, com- 
bined with a specific strain-dependence of the transi- 
tion rates to ensure the proper occupancies of the two 
bound states. A further explanation involved the com- 
pliance of the filaments and the distribution of binding 
sites on the thin filament in addition to a sub-divided 
power stroke l|Huxlev and Tideswellfl996j) . iDukd (l999) 
has suggested that the flat T 2 curve of a muscle fiber can 
arise despite an instability in the dynamics of a single 
pair of filaments, owing to the symmetry of a sarcomere. 
We investigate this possibility further in this article. 



II. METHODS 

Simulation of the stochastic evolution of the system 
was performed using the Gillespie kinetic Monte Carlo 
algorithm, which works as follows. In each simulation 
step the rates of all possible transitions are calculated. 
The time until the next event is chosen as a random num- 
ber with an exponential distribution and the expectation 
value given by the inverse of the sum of all rates. The 
event itself is chosen randomly with a statistical weight 
proportional to its rate. In the situations with stiff (non- 
compliant) filaments and continuous binding sites the 
transition rates can be factorized into factors that depend 
only on the backbone position (which arc the same for all 
motors in a group) and factors that only depend on the 
binding position of a motor (which does not change with 
time unless that motor undergoes a transition) . This al- 
lowed us to use a very efficient {n log n) algorithm based 
on binary trees. 

We assumed complete mechanical relaxation of the 
system in each step, i.e. the strain of all elastic el- 
ements is equilibrated before the next transition takes 
place. The structure involving sarcomeres, filaments and 
myosin heads was described as a circuit of elements with 
given resting lengths and compliances. The strain of ev- 
ery cross-bridge was calculated, given the constraint of 
fixed total length of the system (isometric conditions), 
or of fixed force acting on the ends of the system (iso- 
tonic conditions). 

The T 2 transients were always determined 2 ms after 
the stretch/release. All other parameter values are sum- 
marized in Table [I] 
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FIG. 1 Reaction scheme for a myosin head. The dashed lines 
represent transitions which are considered sufficiently slow to 
be neglected in the calculation. 



III. MOTORS ACTING BETWEEN A SINGLE PAIR OF 
FILAMENTS 

A. Isometric transient in a swinging lever-arm model 

We model the chemical cycle of a myosin head as de- 
scribed in Refs. l)DukeUl999l l2000j) and shown in Fig.[TJ 
A molecule can exist in either a detached state (D) or 
in one of the attached states, AI (with ADPPi) or A2 
(with ADP). T he state AI is a lso referred to as the weak- 
binding state l)Brennerlll99H) . The chemical transition 
between states AI and A2 is concomittant with a con- 
formational change of the molecule, in which the lever 
arm moves through displacement d. From state A2 the 
head undergoes another conformational change with a 
lever-arm shift 5, associated with ADP release. In the 
ATP-rich physiological environment, this step is quickly 
followed by detachment of the head and we therefore skip 
the transitional state (A3) in our model. The attach- 
ment and detachment rates, which are summarized in 
Table U determine the shape of the force- velocity rela- 
tion, but have little influence on the transient response 
to length steps, which mainly depend on the power-stroke 
displacement d and the elastic constant K of the myosin 
cross-bridge. We assume that the transitions Af^A2 
take place on a faster time scale than the detachment 
and reattachment of heads. Therefore, once a head has 
bound with strain £ in the state AI, its state can be de- 
scribed as a statistical ensemble of the states AI and A2 
with probabilities given by the Boltzmann factors 
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where AG stro ke is the free energy change associated with 
phosphate release. 

Now consider a filament containing N myosin motors 
interacting with a actin filament which is held to prevent 
it from sliding. If Na heads are bound to the actin with 
a distribution of strains $(£) (normalized to I), then the 
isometric force that they generate is 



T = N A J + i*(Od)d£ 



(2) 



If a stretch Aa; is suddenly applied to the pair of fila- 
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TABLE I Model parameters 
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Actin elasticity 


7 
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a The numbers in brackets show an alternative parameter set, cor- 
responding to the oscillating case 



ments, the force will first change to 

Ti(Ax) =T + N A KAx , (3) 

as the cross-bridges are deformed. In the next instance, 
transitions between states Al and A2 occur and the prob- 
ability distribution re-equilibrates to P2(£ + Ax). As a 
result the force adjusts to 

T 2 (Ax)=N A J<P{OK{t + Ax + P 2 (Z + Ax)d)dt. (4) 

As a first approximation, we can neglect the distribu- 
tion of strains and assume £ = for all myosin heads, 
which allows us to calculate the T2 curves analytically. 
Their shape depends on two dimensionless parameters, 

Kd 2 

61 = oXr ' (6) 

ZIACr stroke 

which provide a measure of the energy stored in the 
elastic element when the head performs a power stroke 
(Al— >A2). The parmeter eo measures this energy rela- 
tive to the thermal energy, while e\ compares it to the 
magnitude of the chemical free-energy change AG s troke 
that accompanies the power-stroke transition. Note that 
if ei < 1, the power stroke can occur immediately af- 
ter the myosin head binds to the thin filament; but 
if ei > 1, the conformational change is energetically 
inhibited initially, and it is only likely to occur once 



the thin filament has been moved forward by the ac- 
tion of other motors. On grounds of efficiency of en- 
ergy transduction, we expect the value of t\ to be as 
high as possible, while ensuring that there is sufficient 
chemical energy to drive the power stroke, i.e. ti « 1 
llDukeUl999|) . The power-stroke displacement d has been 
measured directly in a number of single-molecule exper- 
iments, and the val ues obtained range betwee n 5 and 
10 nm (reviewed bv iTvska and Warshawl l2002|) . Data 
on the crossbridge elasticity K are less reliable, because 
there a re inevitably other sources of compliance in the 
system l|Veigel et allll998l) . But indirect evidence is pro- 
vided by the en ergetic e f ficiency of muscle, which peaks 
at about 50 % llBarclavt Il998t iKushmerick and Daviesl 
Il969l IPiazzesi et all 120 02 £1 . Thus a lower estimate for 
the energy stored in the elastic element following the 
power stroke is one half of the free-energy change accom- 
panying the hydrolysis of an ATP molecule AGatp ~ 
20keT, which leads to eo ^ 10- Taking into account that 
not all energy stored in the spring can be converted to 
mechanical w ork one obtains a better estimate eo ~ 20 
iDukeL 119991) (which is consistent with ei w 1 if most 
of the energy of hydolysis is used to power the stroke 
AG s trokc ~ AGatp)- This corresponds to an elastic 
constant K — 2.5pN/nm if a power-stroke distance of 
d = 8 nm is assumed. 

The form of the T 2 curve depends critically on the 
values of these dimensionless parameters. 

If eo > 2, an i nterva l with a negative slope (a hystere- 
sis) exists IjHill Il974f) . This means that under isotonic 
(constant load) conditions the system can be bistable. 
The location of the interval of negative slope depends on 
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FIG. 2 The T2 transient response to a step-wise stretch for 
a group of 300 motors in the isometric state for two differ- 
ent values of the parameter ei. The dotted and the dashed- 
dotted lines shows the response of a group of motors which 
are pulling against an elastic element instead of being com- 
pletely isometric before the length step (although the length 
step is then imposed on the motors alone, without the elas- 
tic element). Other parameters were: AGatrokc/ftsT = 15, 
e = 19.3, e 2 = Kd8/k B T = 1.6 (when ei = 1.33) or 
e = 11-6, e 2 = Kd6/k B T = 0.95 (when ei = 0.8), 71 = 2 
and 72 = 0.05. 



the value of t\. If the hysteresis spans the origin, then the 
state with Ax = will be unstable under isotonic con- 
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the freedom to move. This situation occurs, for exam- 
ple, if both filaments are held at their ends by a flexible 
spring (in which case the motors cause the filaments to 
slide until the stall force is reached), or if the filaments 
suffer a constant load which is precisely adjusted to pre- 
vent net sliding. In this case the conditions are actually 
isotonic (or nearly so) and the instability discussed in the 
previous section implies that there is no steady state of 
the system with zero velocity when Eq. (J7J) is satisfied. 
Instead, as shown in Fig. [31 the pair of filaments oscil- 
lates in near-isometric conditions. The mechanism that 
generates the oscillations is the following. An ensemble of 
bound motors pulling against a constant load has two sta- 
ble configurations: one with a majority of motors in state 
A2, the other with the majority in state Al. Because the 
detachment of motors is faster from state A2 than from 
state Al, the total number of bound motors decreases 
when the system is in the first configuration. As a result, 
the load per motor increases until the remaining motors 
cannot support the load any more, whereupon the sys- 
tem flips into the second configuration. Subsequently, the 
number of bound motors starts to grow again, and when 
the load per motor falls below a critical level the system 
flips back to the first configuration. Repetition of the cy- 
cle gives rise to an oscillation whose asymmetry reflects 
the differences in attachment and detachment rates. 

Could the isometric transient be modified if it is mea- 
sured in near-isometric conditions, instead of exactly 
isometric conditions? Fig. ~~| shows the average T2 re- 
sponse in the situation where the filaments are held at 
the ends by an external spring. In the non-oscillating 
case (ex = 1.33, dashed line) the difference between this 
T2 curve and the one in the strictly isometric state is only 
quantitative. The most important result is that the range 
of negative slope in T2 remains. On the other hand, in the 
oscillating case (ei = 0.8, dot-dashed line) the negative 
slope is flattened out because of the broader distribution 
of strains on cross-bridges. 



With a value of e% » 1 the T 2 curve always has a neg- 
ative slope around the stationary point. For example, 
with e = 20, Eq. Q yields 0.85 < e x < 1.22. However, 
it should be bourne in mind that this calculation does 
not take into account the distribution of strains on cross- 
bridges within the ensemble, and so the actual range 
might deviate slightly from this estimate. 

An example of the T2 transient for a group of myosin 
heads between two firmly clamped filaments is shown in 
Fig.H 



B. Group of stiffly coupled motors under near-isometric 
conditions 

For a single pair of filaments, true isometric conditions 
imply that both filaments are held fixed. However, ex- 
perimental conditions are usually near-isometric, in that 
the average velocity is zero, but the filaments still have 



C. Effect of filament elasticity 

The compliance of the thick and thin filaments is also 
expected to affect the T 2 curves. We model it by intro- 
ducing a linear elasticity in the backbone connecting the 
heads (which is, to first order, equivalent t o an elasticity 
in th e track the motors are running on) {Vilfa n et all 
1998). With 7 being the linear modulus of the fila- 
ments, which has t he measured value 7 = 44000 pN 
l|Koiima et HI 0*994), 

the spring constant of a filament 
segm ent of length L is Kt — L~ lr ). In (Vilfa n~et all 
1998) a linear two-state model was used to study the 
effect of filament elasticity on force-velocity relations, 
showing that they reduce the isometric force if the com- 
pliance of a filament segment between two bound mo- 
tors becomes comparable to the motor spring constant 
K . Thin filaments are just stiff enough to prevent sig- 
nificant losses due to this effect. However, the filament 
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FIG. 3 A) A group of motors pulling against a spring. The curve shows the spring extension y as a function of time t. When 
the motors reach stall the average velocity is zero, but fast oscillations occur on the length scale of the power-stroke. The 
observed oscillations are characteristic of the case where the criterion in Eq. £J i s fulfilled. The simulation involved 300 motors 
pulling on an elastic element with spring constant K cxt = lOpN/nm, with other parameters equal those of the dashed line in 
Fig- El B,C) The mechanism of oscillations: (1) When the majority of bound motors are in state A2, the average detachment 
rate is higher than the attachment rate, which leads to a fall in the number of attached motors (2). The force per bound 
motor, shown in (B), therefore increases. When the upper limit of the hysteresis is reached, the system abruptly jumps to the 
other fixed point in which most bound motors are in the state Al (3). Because the detachment rate is low in the state Al, 
the number of bound motors increases again (4) whereby the average force per bound motor drops. This eventually leads to a 
transition back to the fixed point in which most motors are in state A2 (1), at which point the oscillation cycle repeats. 



compliance can have a significantly bigger effect on the 
transient response. The result of a simulation is shown in 
Fig. Q] Although the filament compliance almost halves 
the total stiffness of the system, it leaves the major part 
of the interval of negative slope unaffected. 



D. Effect of discrete binding sites 

So far we have assumed that a myosin head can bind 
anywhere on the actin filament with equal probabil- 
ity. However, in reality the binding sites on actin are 
c = 5.5 nm apart. In addition, the actin filament has 
a helical structure with a half-pitch of about 38 nm (7 
monomers). We therefore propose that the binding rate 
to a particular site is proportional to the Boltzmann fac- 
tor determined by the spring distortion energy, composed 



of a longitudinal and a lateral component 
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with £ — x — (i + 7j)c and (f> — n^, where j denotes 
the repeat on the helix and i the consecutive number of 
the actin monomer. The total binding rate at a given 
position x then reads 
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These curves have b een measured experi mentally using 
SI myosin heads by ISteffen et alJ (§2001) and the data 
were fitted with Ka = 15pN/nm, the value which we 
use here. The position-dependence of the binding rate 
for these parameters is shown in Fig. [SJ A similar dis- 
tribu tion has also measured for myosin V <|Veigel et all 
1200211 . 
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FIG. 4 T2 curves for a stiff (solid line) and an elastic (dashed 
line) filament. The simulations were performed with a group 
of N = 150 myosin heads, attached to a filament of a total 
length 1000 nm and elasticity 7 = 44000 pN. The dotted line 
was obtained when an additional segment of length 500 nm of 
free actin between the clamp and the active part was taken 
into account. The dot-dashed line shows the same situation, 
but additionally taking into account discrete binding sites on 
actin. In the latter case the binding rate fcbind = 100 s -1 
was adjusted to give a realistic fraction of attached heads in 
the isometric state, therefore 71 = 0.8 and 72 = 0.125. All 
these effects together are still not able to completely cancel 
the negative slope. 




FIG. 5 The binding rate fcbind(a:)/fcbind °f a myosin head as 
a function of its position along the actin filament. The longi- 
tudinal spring constant was K = 2.5pN/nm and the lateral 
spring constant was Ka = 15 pNnm. 



The dot-dashed line in Fig. 0] shows the T2 re- 
sponse when taking into account discrete binding sites 
and the elasticity of actin filaments. As noted by 
iHuxlev and Tideswelll l)l996|) . both the filament compli- 
ance and the discrete binding sites contribute to the flat- 
tening of the T 2 curve. However, our simulation indicates 
that their combined effect is not sufficient to cancel the 
negative slope with the parameters used here. We there- 
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FIG. 6 A) A single actin filament interacting with a single 
myosin filament. B) A pair of actin filaments, interacting 
with a myosin filament. C) n = 3 pairs of actin filaments, 
connected in series. D) A muscle-like structure with n — 3 
sarcomeres connected in series and each sarcomere containing 
p = 3 thick filaments, connected in parallel via the Z-discs. 



fore conclude that the absence of negative slope in the 
measured T2 curves cannot be explained in terms of a 
single-filament model, but requires taking into account 
the action of serially connected sarcomeres. 



IV. MOTORS IN MUSCLE SARCOMERES 
A. Redistribution of filaments 

The situation changes essentially if we take into ac- 
count that the experiments are performed on whole mus- 
cle fibers, where several hundreds of sarcomeres are con- 
nected in series, as shown in Fig. The tension in each 
half-sarcomere must be the same, but the extension of 
each half-sarcomere can differ, provided that their sum 
equals the total stretch of the muscle fiber. To see what 
happens in this situation, it is instructive first to look at 
a minimum unit consisting of two filament pairs, joined 
back-to-back as shown in Fig. [SJs. This unit represents 
two thin filaments and a thick filament, the basic build- 
ing block of a sarcomere. A sudden stretch of the unit 
will first displace both halves equally, according to the 
elasticity of cross-bridges. But in the next phase, when 
the distribution of states Al and A2 equilibrates, the 
position of the thick filament can become unstable. This 
means that a slight fluctuation in one direction will cause 
the motors pulling the thick filament in that direction to 
increase their force and those on the other side to de- 
crease it, and the thick filament will jump sideways to 
one of the stable points. A schematic example is shown 
in Fig. □ 
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FIG. 7 The force-displacement relations of the left and the 
right half of a myosin filament before (A) and after (B) a 
length step. The sum of their displacements must equal the 
total sarcomere stretch, Axi + Ax2 = 2Ax. This condition 
is imposed by plotting the force-displacement-relation of one 
set of motors as a function of Axi (a;-axis) and that of the 
other set as a function of Ax2 = 2Ax — Axi (upper a:-axis). 
The forces on the myosin filament are in equilibrium if the 
forces produced on both sides (T2(A:ri) and T2(Ax2)) are 
equal, which is given by the intersection of both curves. The 
equilibrium is stable if the first curve crosses the second from 
below and unstable if it crosses it from above. If there are 
three stationary points (B), the central one is always unsta- 
ble and those to either side stable. The force-displacement 
relation of the two-filament unit is given by the force at the 
stable intersection as a function of the mean displacement Ax. 
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FIG. 8 The T2 response of a single actin-myosin filament pair, 
compared with 2, 4, 8, and 16 filament pairs in series. Each 
myosin filament contains 300 heads. 



The result of a simulation using different numbers of fil- 
ament pairs, connected in series (the situation illustrated 
in Fig. |U) , is shown in Fig. |S| The interval of negative 
slope is increasingly flattened out as the number of pairs 
increases. Apart from the cancellation of the hysteresis, a 
shift in the curves is also observable, corresponding to the 
transition between the strictly isometric conditions and 
those where a filament is held under isotonic conditions 
before the stretch, cf. Fig. [21 The effect of redistribution 
can best be seen if one plots the distribution of stretches 
of individual half-sarcomeres against the total stretch per 
half-sarcomere (Fig.^J. For large stretch amplitudes the 
distribution is concentrated on the diagonal, in agree- 
ment with the expectation that all filaments experience 
the same stretch. For two filament pairs in series, the dis- 
tribution splits in the intermediate range. The two peaks 
represent two stretches for which the single-filament T2 
curves produce the same force. Their position is deter- 
mined by the condition that the average stretch must 
equal the value on the x-axis. With four filament pairs 
in series, a similar situation occurs, but with three differ- 
ent regimes. First three out of four pairs get displaced in 
one direction while the fourth is displaced in the other di- 
rection. Then this ratio changes to two against two and 
finally to one against three, before the regime with all 
filaments being subject to the same displacement settles 
in. With higher numbers of filaments in series (Fig. |5J; 
shows the situation for 32) this pattern becomes increas- 
ingly continuous. In the regime where the negative slope 
stretches over the isometric point, the individual fila- 
ments perform small oscillations (Fig.[3J) while the whole 
muscle is held at a constant length. These oscillations 
blur the distribution further (Fig. . 

T2 curves for four different values of K are shown in 
Fig.UHl 
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FIG. 9 Plot showing the probability density of transient 
stretches of individual half-sarcomeres (y-axis) against the to- 
tal stretch per half-sarcomere. The three plots were obtained 
for (A) 2, (B) 4 and (C) 32 filament pairs in series, each con- 
taining N = 300 motors. (D) shows a simulation with 32 
filament pairs and t\ = 0.8, in which spontaneous oscillations 
of individual filaments take place, thereby blurring the indi- 
vidual shifts. 



B. Isometric state 

It is instructive to take a closer look at the parameters 
of the globally isometric state (by globally we mean that 
the total length of a fiber composed of many sarcom- 
eres is constant, rather than the positions of individual 
filaments) in the two cases: one where the individual fil- 
aments oscillate and one where their positions just fluc- 
tuate. All data were determined in a situation where 
a large number (32 in the simulation) of filament pairs 
were connected in series and the total length held con- 
stant (situation (c) in Fig.^J. One expects the oscillating 
filaments to consume more ATP than the non-oscillating 
filaments. The results are summarized in Table ITT1 An 
exact comparison with experimental data is difficult, be- 
cause the fraction of attached heads in the isometric state 
is not known reliably, and also because the experiments 
were carried out at different temperatures and with dif- 
fer ent types of myosin. F o r exam ple, the value measured 
bv lBarsotti and Ferenczil |l988). 0.4 ATP molecules per 
second per head (with 25% of the heads attached this 
gives 1.6 s _1 per attached head), is compatible with the 
first sce nario. However, other m easurements give higher 
values l)Ebus and Stienenl . Il996ft . more consistent with 
the second case. We therefore conclude that from the 
experimental values of the ATPase rate in the isometric 
state neither scenario can be excluded. 

At high values of ei, the model shows an interesting 



K=1.5 pN/nm 
K=2.0 pN/nm 
K=2.5 pN/nm 
K=3.0 pN/nm 




FIG. 10 T2 curves with filament redistribution. The simu- 
lations were performed with n = 16 filament pairs in series, 
each containing N = 300 myosin motors. Different curves 
show data for different values of the cross-bridge elasticity 
K, corresponding to values of ei between 0.8 and 1.6, and eo 
between 11.6 and 23.2. 



TABLE II The parameters of the globally isometric state in 
the non-oscillating and the oscillating case. 



Elastic constant 

eo 

ei 

02 



Non-oscillating Oscillating 
K = 2.5pN/nm K = 1.5pN/nm 
19.3 11.6 
1.33 0.8 
1.6 0.95 



Fraction heads attached 

! I ' ■: nls in state A2 
Attached heads 
Force 



attached head 
ATPase 

attached head 



92% 
3.5% 
4.8 pN 
1.7s- 1 



68% 
32% 
7.4 pN 
17.4s -1 



feature. In the isometric state the number of cross- 
bridges in the state A2 can be very low, although the 
generated force per attached cross-bridge is as high as 
5 pN, in agreement with single-filament measurements 
ijKawai et all 12000). This allows the muscle to support 
a load un der isometric conditions with little ATP con- 
sumption llDukei . 120001) . It might at first sound para- 
doxical that most of the isometric force is generated by 
the pre-power-stroke state Al. This is possible because 
the myosin heads bind with a stochastic distribution of 
strains. Because those with a negative strain are more 
likely to undergo the power-stroke and then detach, the 
remaining ensemble produces a positive force. In this as- 
pect, the power-stroke actually serves to eliminate nega- 
tively strained cross-bridges rather than direct generation 
of positive force. This notion is in agreement with recent 
experiments, which have shown that high phosphate con- 
centration does reduce the isometric force significantly 
ijCooke and Patej. Il985() , but it does not have a visible 
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effect on the conformati on of the myosin h eads or even 
their catalytic domains ijBaker et all Il999|) . But let us 
stress again that this holds only for the isometric state. 
The force in a contracting muscle originates mainly from 
the state A2. 



C. Isotonic response 

Another important class of experiments which provides 
information on the actomyosin interation involves the iso- 
tonic transient. Here the applied force is initially set 
at the value of the stalling load T , so that the fiber 
is prevented from contracting. The force is then sud- 
denly changed and held constant at a different value, 
while the length of the fiber is recorded. Some early ex- 
periments showed that following a small step change of 
load, damped oscillations were impose d on the steady 
contraction or e xtens ion of the fiber ijGranzier et all 
1199ft iPodolskvl ^60). Such oscillations are particu- 
larly clear in recent exper iments on single muscle fibers 
l|Edman and Curtinl EkH^ . 

One possible cause of this oscillatory response has pre- 
viously been suggested on the basis of the stochastic 
model of the actomyosin interaction used in this arti- 
cle. When ei > 1, the chemical cycles of myosin motors 
on the same filament can become synchroniz ed at loads 
close to the stalling force (Duke, 1999, 2000). A pair of 
filaments then slides in a step- wise fashion under isotonic 
conditions. But during steady shortening, the motors on 
different filaments within the same muscle fiber operate 
out of phase, so that there is no macroscopic manifesta- 
tion of the steps. However, an abrupt change in the load 
can cause the synchronization of a large fraction of the 
bound motors, w hereupon the steps do become observ- 
able l|Dukel . ll999|) . Because the correlation of the motors 
soon decays, the macroscopic steps fade and a damped 
oscillation is seen. 

A much stronger oscillatory response is seen in the 
regime where individual pairs of filaments perform os- 
cillations in the near-isometric state, as described in 
Sect. 1111. Bl The synchronized oscillations can then be 
very pronounced following a small decrease in the load, 
as shown in Fig.lllb. On the other hand, no damped os- 
cillations are observed after a larger drop in the load, e.g. 
to To /3 (Fig. II lb ) , because the individual filament pairs 
immediately move out of the hysteretic regime. These 
properties are in agreement with recent experiments on 
th e isotonic response of sing le muscle fibers carried out 
bv lEdman and Curtinl 1)200i|k 

In steady isotonic conditions, another kind of insta- 
bility can arise due to a hysteresis in the force-velocity 
relationship. The possibility of such an instability was 
first discussed in the conte xt of a two-state ratchet model 
i|Jiilicher and ProstL Il995|) and subsequently in a kinetic 
cross -bridge mode l with a strain-dependent detachment 
rate (jVilfan et all Il999(l . With the parameters we use 
this hysteresis, though existent, covers a rather small 




FIG. 11 The isotonic transient response of a sarcomeric struc- 
ture to: (A) a small drop, (B) a large drop, and (C) a rise in 
the load, when the initial value of the load is chosen precisely 
to stall the contraction. The upper graph shows the force 
per myosin filament as a function of time. The lower graph 
displays the length change per half-sarcomere and the num- 
ber of attached myosin heads per myosin filament, which is a 
measure of the fiber stiffness. The data were computed with 
N = 300, ei = 0.8 and with 50 myosin filaments acting in par- 
allel. The results were averaged over 125 events (which has 
the same effect as simulating that number of sarcomeres in se- 
ries). A small drop in the load synchronizes the cross-bridges 
and therefore causes observable macroscopic oscillations (A). 
The oscillations are less pronounced, but still visible following 
a small increase in the load (C). 
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range of velocities and onl y results in a sma ll inflexion 
in the force- velocity curve jfDukd. 119991 12000^ . 



V. DISCUSSION 

In the past, the isometric transient model of muscle 
has been modelled by considering the dynamics of a sin- 
gle pair of filaments. A significant problem with this 
approach is that, with a general choice of param eters 
in a cross-bridge model, the T2 curve is not flat l)Hill 
Il974|) — it typically has either a positive or a negative 
slope for limitingly small step displacements. Thus, in 
order to reproduce the experimental T2 curve of mus- 
cle, some of the paramet e rs ha v e needed to be finely 
adjus ted ijEisenberg et all 1198(1 iHuxlev and Simmons, 
Il97l|) . and limitations have been imposed on the size 
d of the power-stroke and the rigid ity K of the cross- 
bridges l|Huxlev and Tideswell fl996j) / In this article we 
have investigated the isometric transient response using 
a stochastic model of the actomyosin cycle. When values 
of the parameters d and K are chosen to explain other 
characteristics of muscle, like the force-velocity relation 
and the efficiency of energy transduction, the Ti curve 
of a single pair of filaments displays a region of negative 
slope, which cannot be eliminated by factors including 
filament compliance and the discrete nature of the bind- 
ing sites on the actin filament. However, we argue that 
the symmetric structure of sarcomeres must be taken into 
consideration when computing the isometric transient re- 
sponse of muscle. Following a step change of length of a 
muscle fiber, a redistribution of sarcomere lengths occurs 
within the fiber. Some half-sarcomeres contract, while 
others extend. This redistribution always eliminates the 
negative slope, leading to a flat Ti curve. It is a generic 
feature of unstable elements connected in series, and does 
not require any special values of model parameters. 

Our model shows that there are two different regimes 
of the microscopic dynamics in near-isometric conditions. 
For a range of values of the parameter ei close to unity, 
the isometric point falls in the interval where the slope of 
the Ti curve is negative. In this case individual filament 
pairs oscillate with small amplitude. In the other regime, 
where the T 2 curve has positive slope at the isomet- 
ric point, the individual filaments are stationary, apart 
from stochastic fluctuations. The macroscopic manifesta- 
tions of these two regimes differ in few respects. Because 
the oscillations of different filament pairs have different 
phases, oscillatory motion is not normally observable on 
the scale of a whole muscle fiber in steady conditions. 
However, a sudden change of load can synchronize the os- 
cillations and thereby make them visible. The existence 
of damped oscillations in the isotonic transient re sponse 
of single muscle fibers (|Edman and Curtini 120011) there- 
fore argues in favor of the oscillating regime. We note, 
however, that damped oscillatio ns can also b e a manifes- 
tation of step-wise shortening llDukel. QUI, which can 
exist in both regimes. Because the efficient transduction 



of energy demands t\ ~ 1, which is close to value of 
this parameter at the boundary of the two regimes, it is 
possible that both regimes exist depending on conditions 
such as the myosin isoform, phosphate concentration, pH, 
io nic strength and temperatu re. Indeed, measurements 
bv lEdman and Curtini (|200l|) show a dependence of the 
oscillation decay on the solution pH and on muscle fa- 
tigue. Further experiments, in which conditions are sys- 
tematically varied, could shed more light on the mecha- 
nism of oscillation. 

As a final remark, we emphasize that according to our 
model, the T2 curve of an individual pair of filaments 
differs from that of a muscle fiber. Recently, assays have 
been developed to measure the force-veloci ty relation of 
a sin gle filament within a half-sarcomere ijKawai et all 
2000). In order to test our predictions, it would worth- 
while to develop such high-precision techniques to mea- 
sure the transient response of a single filament. 
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